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Abstract 

We study large deviations in the Langevin dynamics, with damping of order e~^ and 
noise of order 1, as e f 0. The damping coefficient is assumed to be state dependent. 
We proceed first with a change of time and then, we use a weak convergence approach 
to large deviations and their equivalent formulation in terms of the Laplace principle, 
to determine the good action functional. 

Some applications of these results to the exit problem from a domain and to the 
wave front propagation for a suitable class of reaction diffusion equations are considered. 


1 Introduction 


For every e > 0, let us consider the Langevin equation 

r(t) = Kq^{t)) - LMhlgqp + a{q^{t))B{t), 
< ^ 

q^{0) = q e gPO) = p G M'^. 


( 1 . 1 ) 


Here B{t) is a r-dimensional standard Wiener process, defined on some complete stochastic 
basis (U, A", {Jq}) IP)- In what follows, we shall assume that b is Lipschitz continuous and 
a and a are bounded and continuously differentiable, with bounded derivative. Moreover, 
a is invertible and there exist two constants 0 < ao < ai such that oq < ot{q) < ai, for all 
g G Equation 0 can be rewritten as the following system in 

' q^it)=p%t), gqO)=gGM^ 

pgt) = b{q^{t)) - LMLI+ a{q^{t))B{t), f{0) = p G M", 

^ 6 
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and, due to our assumptions on the coefficients, for any e > 0, T > 0 and A: > 1, the system 
above admits a unique solution z'^ = G C([0, T]; which is a Markov 

process. 

Now, if we do a change of time and define q^it) := q^{t/e), f > 0, we have 
e^qe{t) = b{qe{t)) - a{q^{t))q^{t) + y/ea{qe{t))w{t), 

V 

q,{0) = g E M"*, g,(0) = - E M"*, 

V 6 

where w{t) = f > 0, is another M^-valued Wiener process, defined on the same 

stochastic basis (fl, J", P). 


In the present paper, we are interested in studying the large deviation principle for 
equation (1.2), as e | 0. Namely, we want to prove that the family {ge}e>o satisfies a 
large deviation principle in C{[0,T]', H), with the same action functional I and the same 
normalizing factor e that describe the large deviation principle for the first order equation 


9e{t) 


Kaejt)) 

a{9e{t)) 


+ Ve 


^{9e{t)) 

0^{9e{t)) 


w{t), 


g,{0) = g E M'^. 


(1.3) 


In particular, as shown in Section this implies that the asymptotic behavior of the 
exit time from a basin of attraction for the over damped Langevin dynamics( |l.l[ ) can be 
described by the quasi potential V associated with /, as well as the asymptotic behavior of 
the solutions of the degenerate parabolic and elliptic problems associated with the Langevin 
dynamics. 

Moreover, in Section we will show how these results allow to prove that in reaction- 
diffusion equations with non-linearities of KPP type, where the transport is described by 
the Langevin dynamics itself, the interface separating the areas where is close to 1 and 
to 0, as e 4 , 0, is given in terms of the action functional I, as in the classical case, when the 
vanishing mass approximation is considered. 


In [8] and [3], the system 


w'tlA*) = Kqi^A^)) - Aq^lA^^))Q^lAt) + A^c^iq^lAt))w{t), 

q^A^) = q£R^, = - E M'^, 

V e 


for 0 < ;U, e << 1, has been studied, under the crucial assumption that the friction coefficient 
a is independent of q. 

It has been proven that, in this case, the so-called Kramers-Smoluchowski approxi¬ 
mation holds, that is for any fixed e > 0 the solution of system (1.4) converges in 
(^([O,r ];as fj, 0, to g^:, the solution of the first order equation (1.3). More¬ 
over, it has been proven that, if VAq,p) is the quasi-potential associated with the family 
{Q'/i,e}e> 0 ) for /i > 0 fixed, then 


lini inf VAq,p) = y{q), 

11^0 p 6 Rd 
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where V is the quasi-potential associated with the action-functional I. 


In [^, equation (1.4) with non constant friction a has been considered and it has been 
shown that in this case the situation is considerably more delicate. Actually, the limit of 
qfi,e to has only been proven via a previous regularization of the noise, which has led 
to the convergence of to the solution of the first order equation with Stratonovich 
integral. 

Finally, we would like to mention that in the recent paper by Lyv and Roberts, an 
analogous problem has been studied for the stochastic damped wave equation in a bounded 
regular domain D C with d = 1,2, 3, 

d‘^u(t,x) ^ , ,, du(t,x) ^dw(t,x) 

e- =/2.u{t,x) + f{u{t,x)) -^ ' 


dt'^ 


dt 


dt 


u{t,x) = 0, x^dD, tt(0, x) = tto(x), — ^ — '^-=vo,{x) 


where e > 0 is a small parameter, the friction coefficient is constant {a = 1), w{t,x) is 
a smooth cylindrical Wiener process and / is a cubic non-linearity. By using the weak 
convergence approach, the authors show that the family {ne}e>o satisfies a large deviation 
principle in C{]fi,T\;L?‘{D)), with normalizing factor and the same action functional 
that describes the large deviation principle for the stochastic parabolic equation. 


As mentioned above, in the present paper we are dealing with the case of non-constant 
friction a and fi = e^. Dealing with a non-constant friction coefficient turns out to be 
important in applications, as it allows to describes new effects in reaction-diffusion equations 
and exit problems (see section |^. Here, we will study the large deviation principle for 
equation( |1.2[ ) by using the approach of weak convergence (see [T] and [2]) and we will show 
the validity of the Laplace principle, which, together with the compactness of level sets, is 
equivalent to the large deviation principle. 

At this point, it is worth mentioning that one major difficulty here is handling the 
integral 

/ exp ( - / a{qe{r))dr] a{q^{s))dw{s), 


and proving that it converges to zero, as e j, 0, in L^{Q; (^([O, T]; Actually, as a is non¬ 
constant, the integral above cannot be interpreted as an Ito’s integral and in our estimates 
we cannot use Ito’s isometry. Nevertheless, due to the regularity of (?e(t), we can consider 
the integral above as a pathwise integral, and with appropriate integrations by parts, we 
can get the estimates required to prove the Laplace principle. 


2 The problem and the method 

We are dealing here with the equation 

e^qe{t) = b{qe{t)) - a{q^{t))q^{t) + y/ea{qe{t))w{t), 

V 

q,{0) = qeR‘^, qe{0) = - £ R^, 

V 6 
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Here w{t), t > 0, is a r-dimensional Brownian motion and the coefficients b, a and a satisfy 
the following conditions. 

Hypothesis 1. 1. The mapping b \ ^ is Lipschitz-continuous and the map¬ 
ping cr : —)• is continuously differentiable and bounded, together with its 

derivative. Moreover, the matrix (7{q) is invertible, for any q G and —)• 

£(]R'’,]R'^) is bounded. 

2. The mapping a : —>• M belongs to (M'^) and 


inf a{x) =: oq > 0. 


xei 


( 2 . 2 ) 


In view of the conditions on the coefficients a, b and a assumed in Hypothesis for 
every fixed e > 0, equation (2.6) admits a unique solution = {qe,Pe) £ £^(0,T;M‘^), with 
r > 0 and A: > 1. 

Now, for any predictable process u taking values in L^([0,T];M^), we introduce the 
problem 

rr(n^(+\\ 

(2.3) 


5“(0) = g G M'". 


a{g^{t)) a{g^{t)) 

The existence and uniqueness of a pathwise solution to problem (2.3) in (^([O,T];M'^) 
is an immediate consequence of the conditions on the coefficients b, a and o. that we have 
assumed in Hypothesis [Tj 

In what follows, we shall denote by Q the mapping 

Q : L2([o,r];M’') ^ C([0, T]; u ^ Q{u) = g^. 

Moreover, for any / G C([0,T];M'^) we shall define 


^(/) = 


\u{t)\‘^dt : f = Q{u), uG L2([0,r];M^) 


with the usual convention inf 0 = +oo. This means that 


Iif) = 


1 


o^if{s))cr (/(s)) /(s) - 


Kfjs)) 

a{f{s)) 


ds, 


(2.4) 


for all / G IH^’2(o,r;M'^). 

If we denote by g,, the solution of the stochastic equation 

= + 9,(0) 

a{ge[t)) a{ge{t)) 


(2.5) 


we have that I is the large deviation action functional for the family {ge}e>o in the space 
of continuous trajectories C([0,T];M'^) (for a proof see e.g. [H]). This means that the 
level sets {/(/) < c} are compact in (£([0,T];M'^), for any c > 0, and for any closed subset 
F C (£([0, T]; M*^) and any open set G C (^([O, T]; M*^) it holds 


limsupelogP( 5 e G F) < —I{F), 
£-(•0 + 


liminf elogP(ge G G) > —1(G), 

e^0+ 
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where, for any subset A C (^([O,T];M'^), we have denoted 

The main result of the present paper is to prove that in fact the family of solutions 


of equation (1.2) satisfies a large deviation principle with the same action functional I that 


describes the large deviation principle for the family of solutions of equation (2.5). And, 


due to the fact that q^{t) = qe{et), t > 0, this allows to describe the behavior of the over 
damped Langevin dynamics 0 (see Section]^ for all details). 

Theorem 2.1. Under Hypothesis^ the family of probability measures {C{qe)}e>o, in the 
space of continuous paths (^([O, T]; 1^), satisfies a large deviation principle with action func¬ 
tional I. 


In order to prove Theorem 2.1 we follow the weak convergence approach, as developed 
in [T], (see also [2]). To this purpose, we need to introduce some notations. We denote by 
Vt the set of predictable processes in L‘^{Q x [0,T];M''), and for any T > 0 and 7 > 0, we 
define the sets 


5^=<j/G L^(0,r;: 




\f{s)\‘^ds < 7 


= {u £ Vt '■ u £ P — a.s.} . 

Next, for any predictable process u taking values in L^([0, T]; M^'), we denote by qf{t) 
the solution of the problem 

' =biQeit)) - a{qf{t))q^{t) + y/ea{qf{t))w{t) + a{qf{t))u{t), 

( 2 . 6 ) 

g“(0) = g G M'^, g“(0) = ^ G 

V 6 

As well known, for any fixed e > 0 and for any T > 0 and k > 1, this equation admits a 
unique solution qf in L^(n; (^([O, T]; M'^)). 

By proceeding as in the proof of [H Theorem 4.3], the following result can be proven. 

Theorem 2.2. Let {ue}e>o be a family of processes in Sj, that converge in distribution, as 
e i 0, to some u £ S^, as random variables taking values in the space L‘^{0,T;W^), endowed 
with the weak topology. 

If the sequence {g“'}e>o converges in distribution to g^, as e 0, in the space of contin¬ 
uous paths (^([O,T];M'^), then the family {£(ge)}£>o satisfies a large deviation principle in 
(^([O, Tj; M*^), with action functional I. 

Actually, as shown in [2], the convergence of qf^ to implies the validity of the 
Laplace principle with rate functional I. This means that, for any continuous mapping 
A : C([0,r];M^) it holds 

lim —elogE exp f — - A(ge) ) = inf (A(/) + /(/)). 

^ V e ^ V /eC([0,T];R'i) 

And, as the level sets of I are compact, this is equivalent to say that {£(ge)}e>o satisfies a 
large deviation principle in (^([O, T]; M'^), with action functional I. 
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3 Proof of Theorem 2.1 


As we have seen in the previous section, in order to prove Theorem 2.1, we have to show 


that if {ne}e>o is a family of processes in that converge in distribution, as e 0 , to 
some u G 5^, as random variables taking values in the space L^( 0 ,T;M'^), endowed with 
the weak topology, then the sequence {(/“''}e>o converges in distribution to 5 “, as e | 0 , in 
the space (^([O,T];M'^). 

In view of the Skorohod representation theorem, we can rephrase such a condition in 
the following way. On some probability space {Cl, !-,¥), consider a Brownian motion wt, 
t >0, along with the corresponding natural filtration Moreover, consider a family 

of {.Fjj-predictable processes {ue, ft}e>o in L‘^{Ci x [0,T];]R'^), taking values in S^, IP-a.s., 
such that the joint law of {u^, u, w), under IP, coincides with the joint law of u, w), under 
P, and such that 

P — a.s. (3.1) 


lim Ue = u, 
e^-O 


as L^(0, T;M)-valued random variables, endowed with the weak topology. Let g"' be the 


solution of a problem analogous to ( 2 . 6 ), with u and w replaced respectively by Ue and w. 
Then, we have to prove that 


lirn =g 

e^O 


— a.s. 


in (^([O, r];M“). In fact, we will prove more. Actually, we will show that 

limE sup = 0 . 

[0,T] 


(3.2) 


In order to prove (3.2), we will need some preliminary estimates. For any e > 0, we 
define the process 


if,(t) = ^/ee-^'W f e^^^^^a{q^{s))dw{s), t>0. 

Jo 


(3.3) 


Lemma 3.1. Under Hypothesis^ for any T > 0, /c > 1 and 7 > 0, there exists eo > 0 
sueh that for any u G and e G ( 0 , eo] 


k k 


supE\H,{t)\" < Ck,^iT){\q\^ + \p\'^ + l)e^ + Cfc e^ t2e 

S<t 


Moreover, we have 


E sup \H,:{t)\<y/ec^{T){l + \q\ + \p\). 
te [o,r] 


(3.4) 

(3.5) 


Proof. Equation ( |2.6[ ) can be rewritten as the system 
qe^{t)=pf{t), qf{0)=q 

.- 2-11 / 


^pf{t) = h{qf{t)) - a{qf{f))pf{t) + yfea{qf{t))w{t) + a{qf{t))u{t), pf{0) = 


6 







Thus, if for any 0 < s < f and e > 0 we define 

1 

Ae{t,s):=-^J a{q'^{r))dr, A^{t) := A^{t,0), 

we have ^ 

p^{t)= ^e-^'Wp+4 / e-^‘^^’^h{q^{s))ds 
e e Jo 

+ 4 [ e~^J^’‘'^a{q^{s))ue{s) ds + 

Jo e 

Integrating with respect to t, this yields 

q^(t)= q+^ [\-^J^^pds+\ f r e-^J^’^hiqJirAdrds 

c Jo Jo Jo 

+ 4 f f Me(r) dr ds + 4 [ H^{s)ds. 

^ Jo Jo f Jo 

Thanks to the Young inequality, this implies that for any t G [0, T] 

\qJ:{t)\ <\q\+e\p\+cJ^{l + \qj;{s)\)ds + |Me(s)|ds + 4 ^ |ide(s)|ds 

< c^(r)(|g|+e|p| + 1 ) + 4 / \H^{s)\ds+ f \qj{s)\ds, 

e Jo Jo 

and from the Gronwall lemma we can conclude that 

k“(t)| < c^iT) (1 + \q\ + \p\) + c(r)4 r ms)\ds. 

e Jo 

This implies that for any k > 1 


(3.6) 


(3.7) 


\qnt)r<Ck,^{T){\qf + \pf + l) + c,,,iT)e-^^ / ^3)^ ds, eG (0,1]. (3.8) 

Jo 


Now, due to (3.6), we have 


1 _“oG 


go(*-s) 


|p“'=(t)|<-e ^|p| + ^ / e (l + |g“(s)|)ds 

e ^ Jo 

1 /"^ ao(t-«) 1 

+ ^ / e |u,(s)| ds + ^|ide(t)|, 

£ do ^ 


so that, thanks to (3.8), for any e G (0,1] we get 


|p“(t)| < -e ^jpj+c.y(T)(|g| + jpj + 1) + 4 / e |u,(s)| ds + c(r)4|dd6(t)| 
e do e 


(3.9) 
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As well known, if / G and g G (^([Ojf]), then the Stiltjies integral 


f{s)dg{s), t > 0, 


is well defined and, if ( 7 ( 0 ) = 0, the following integration by parts formula holds 

/ f{s)dg{s)= [ {g{t) - g{s))h'{s)ds + g{t)h{Q), t>0. 

Jo Jo 


(3.10) 


Now, the mapping 


[0,+oo) s ^ 


is differentiable, P-a.s., so that the stochastic integral in (3.3) is in fact a pathwise integral. 
In particular, we can apply formula (3.10), with 


h{s) = g{s) = w{s), 


and we get 


Heit) = Ve iw{t) -wis))e ds 


(3.11) 


+^/iu;(t)e-^^Wc7(g). 


Thanks to (3.9), this yields for any e G (0,1] 






< Cjy 


.{t)\ < c^/e f \w{t) - tc(s)|^^-2-(1 + ds + c^/e\w{t)\e 

Jo ^ 

f{T){W\ + IpI + l)Ve / ' h(t) - w{t - e^s)|e"“o^ ds 
Jo 

r\ 

+\/eCry{T) / |t(;(t) — ri;(t — e^s)|e““°®|//e(t — e^s)| ds + c-v/e |r(;(t)|e 

Jo 

and hence, for any A: > 1, we have 

< Ck,^{T){\q\^ + N" + l)et r" \w{t) - w{t - e^^)] ds 

Jo 


CQt 


+£2 Cfc^.y(T) / |tt;(t) — r(;(t — e^s)|^e “°^|dde(t — e^s)|^ ds + Cfc e2 |t(;(t)|^e 


kaQt 
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By taking the expectation, due to the independence of \w{t) — w{t — e^s)| with — e^s)| 
and /g ^ ^ \ Hf:{r)\^ dr, this implies that for any e G (0,1] 

t 

1 1 1 / 7^ k 

<Ck^^{T){\q\'^ + \p\'^ + l)e^ / s^e-^^^ds 

Jo 

t 

"^k I 72^ k ^ T k k kOLr\t 

+e^ Cfc^g,(T) / s2e“"o^E|ffe(i — ds + t^e~ 

Jo 

< Cfc,g,(r)(|g|^ + IpI^ + l)e^ +Cfcettte ^ Ck,.y{T) supE |iJ,(s)|''. 

S<t 


Therefore, if we pick cq G (0,1] such that 

e^Ck,^{T) < 


we get (3.4). 


Now, let us prove (3.5). From (3.11), we have 


<y/ec sup |tc(t)| 
te [o,r] 


1 + 


Qn(*-2) 


Ip7'(s)I Js 


< 'Jec sup \w{t)\ 
te [o,T] 




\p^^is)\^ds 


and hence 


E sup |Ffe(t)| <-v/ec(r) [ 1 + e [ E / |p“"(s)|^ds 
te [o,T] \ \ Jo 


Thanks to (3.9), as a consequence of the Young inequality, we get 

rt 


/o 


\pf{s)\^ds<c,{T){l + \q\^ + \p\‘^) + \c{T) [ ms)\^ds, (3.12) 


/o 


so that 


E sup |Ff£(t)| < Vecg,(r)(l + 1^1 + IpI) + ^ c(r) r / E\H^{s)f ds 
te[o,T] \Jo 


Therefore, (3.5) follows from 


□ 


Lemma 3.2. Under Hypothesis\^ for any T > 0, k > 1 and 7 > 0 there exists cq > 0 sueh 
that for any u £ and e G (0, cq) we have 


E sup |g“(t)|^<cfc,^(r)(|g|^ + |p|^ + l)e-t+cfc,.,(r)e2-f. (3.13) 

is [0,T] 
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Proof. Estimate (3.13) follows by combining together (3.4) and (3.8). 


□ 


Now, we are ready to prove (3.2), that, in view of Theorem 2.2, implies Theorem 2.1 


Theorem 3.3. Let {ne}e>o be a family of predietable processes in that converge P-a.s., 
as e 4. 0, to some u € S^, with respect to the weak topology o/L^(0, Then, we have 


lim E sup \qf^{t)-g'^{t)\=0. 
te [o,T] 


(3.14) 


Proof. Integrating by parts in (3.7), we obtain 

?."■ (*)=?+ *+s. w. 


/o Oi{qe^{s)) Jo a{qf^{s)) 


where 


e Jo a[qe\t)) Jo Jo a[qP{s)) 




dw{s) 


^lo a^iql^is)) ds 

(t) + / 2 ( n,( (Va(g“'(s)),p“^(s)) ds =: 

Jo “ we (SjJ ^ 


u. /j_\ \ 


a{qP{t)) 

This implies that 

qf^{t)-g'^{t) = 


\b{qp{s)) 

b{g^{s))] 

ds+ f 
Jo 

W{Qp{s)) 


[a(<7/^(s)) 

a{g^{s)) 

[a{qp{s)) 

a{g^{s))_ 


Ue{s) ds 


+ / /'^XN [Ueis) - n(s)] ds + Re{t). 


Jo a{g^is)) 

(3.15) 

Due to the Lipschitz-continuity and the boundedness of the functions a and 1/a, we have 
that a/a is bounded and Lipschitz continuous. Then, as G Sj, we obtain 


< c 


lo a{g^{s)) 


Ue{s) — u(s)] ds 


2 rt 

2 


+ c|i?,(t)|^+c(r) / \qf^is)-g^isrds 
Jo 


+c{T) [ ke'is) - g''{s)\'^ ds ( [ |u,(s)pds+ sup \g^{s)f 
Jo \Jo se [o,t] 


< c 


'■*cj(5“(s)) |DaM2 '■* 

[ne(s) — u(s)J as 


lo Oi{g^{s)) 


+ c\RS)\'^ + c^{T) I \qf^{s)-g'^{s)\^ds. 


lo 
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By the Gronwall lemma, this allows to conclude that 


sup \q-^{t)-g-{t)\ 
tG [o,r] 


< c^{T) sup 
te [o,T] 


crig^s)) 


lo oi{g^{s)) 
Now, for any e > 0, we define 


[ue(s) — n(s)] ds 


+ C^{T) sup |-Re(t)|. 
te [o,T] 


For any 0 < s < t we have 

r.(t) - r,(s) = ^ K(^) - “W] dr, 

so that, as and u are both in S^, 

|re(t) — re(s)| < CjVt — s, e > 0. 


(3.16) 


As re(0) = 0, this implies that the family of continuous functions is {re}e>o is equibounded 
and equicontinuous, so that, by the Ascoli-Arzela theorem, there exists i 0 and v £ 
(^([O, T]; such that 


lini sup -v{t)\ =0, 

t£ [o,r] 


— a.s. 


On the other hand, as (3.1) holds, for any h £ we have 


lim h) = lim ( — u, —r~~ 7 ^h) 


= 0 , 


L2(0,T;R'') 


SO that we can conclude that v 


0 and 

lim E sup |re(t)| = 0. 
•=“^0 te [o,r] 


Thanks to (3.16), this implies that 


limsupE sup |(jr“'(t) — g'“(t)| < c limsup E sup |i?e(t)|, 
e-5-0 te[o,r] £^"0 te[o,T] 


so that (3.14) follows if we show that 


lim E sup |i?e(t)| = 0. 

t£ [0,T] 


(3.17) 


We have 




\P\ 



e-^ds) 



e ds < c\p\ e. 


(3.18) 
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Moreover 








^b{q^^{s))ds 


< c 



ao(t-s) 

^2 


il + \q^^is)\)ds<ce^ 


1+ sup \q^‘{t)\ 

, [0,T] 


Thanks to (3.13), this implies 


E sup |4^(t)| < c^(r)(|p| + |g| + l)e2, eG(0,l]. 

t& [0,T] 


Next 


E sup l/g (t)| = \/iE sup 
te [o,r] te [o,t] 


0 


dw{s) 


< c(T) V~e. 


Concerning we have 


(3.19) 


(3.20) 





^-A^{s,r 


^b{q'^‘{r)) dr 


a^q^^is)) 


(Va(g“'(s)),p“'(s)) ds 


< 6^0 1+ sup |g^(t)|) / \p^,^{s)\ds, 

V ie [0,T] 


so that, due to (3.13) we obtain 


E sup |/(^(t)| < e2c.^(r)(|g| + IpI + l)e ^ Te / |p“'(s)pds 
te [o,T] V Jo 


As a consequence of (3.4) and (3.12), this yields 


E sup |/,^(t)| < ec^(r)(|gp + |p|2 + 1), eG(0,eo]. 

te [o,r] 


Concerning according to (3.5) we have 


E sup |/,®(t)|<cE sup \H^{t)\ < ^/ec^{T){l + \q\ + \p\). 
t£ [o,r] t£ [o,T] 


(3.21) 


(3.22) 


Finally, it remains to estimate We have 

rt I 


\int)\ = 


lo a^qe^is)) 


H,{s) (Va(g^(s)),p“'(s)) ds 


<c \H^{s)M‘{s)\ds, 
Jo 


so that 


E sup |/®(t)| < c 
te [o,r] 


K\H^{s)\‘^ds [ E\p] 

0 Jo 


(s)p ds 


1 

2 
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By using (3.12), this gives 


E sup \I^{t)\<c^{T){l + \q\ + 
te [o,r] 


E\H,{s)f dsj + ^ E\H,{s)fds, 


so that, from (3.4) we get 


E sup \I^{t)\<ec^{T){l + \q\ + \p\), eG(0,eo]. 

te [o,r] 


This, together with (3.18), (3.19), (3.20), (3.21) and (3.22), implies (3.17) and (3.14) follows. 

□ 


4 Some applications and remarks 

Let G be a bounded domain in with a smooth boundary dG. We consider here the exit 
problem for the process q^{t) defined as the solution of equation (1.1). For every e > 0 we 
define 

:= min{f > 0 : q^{t) ^ G}, := min{f > 0 : qe{t) ^ G}, 


where q^{t) = q’^{t/e) is the solution of equation (2.6). It is clear that 

= 1^, q\r^) = 

In what follows, we shall assume that the dynamical system 

q(t) = b{q{t)), t > 0, 

satisfies the following conditions. 


(4.1) 


Hypothesis 2. The point O ^ G is asymptotically stable for the dynamical system (4.1) 
and for any initial condition q £ 

lim q{t) = O. 

t^OO 

Moreover, we have 

{b{q),i^{q)) > 0, q£dG, 
where v{q) is the inward normal vector at q € dG. 

Now, we introduce the quasi-potential associated with the action functional I defined 
in ([2^ 

V{q)=mi {/(/), /G G([0,T];M'='), /(O) = O, f{T)=q, T > o} 


= - inf ' 


a{fis))a ^(/(s)) (f{s) - 


10 


«(/(«)) 


ds, /(O) = O, /(T) = q,T>0}. 


It is easy to check that, under our assumptions on a{q), the quasi-potential V coincides 
with 

^ inf I f cr“^(/(s)) (^f{s) - a{f{s))b{f{s))^ ds, /(O) = O, f(T) = q, T > o| . (4.2) 
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Theorem 4.1. Under Hypotheses^ and[^ for each q £ {q £ G : V{q) < Vq} and p G 
we have 


and 


where 


lmelogE(,_p)r" = Im elogE(g_p)re = Vq, 


lim elog H = lim elog = Vq, in probability, 

e —>0 e —^0 


(4.3) 

(4.4) 


Vb := min y(g). 
qS oG 

Moreover, if the minimum of V on dG is achieved at a unique point q* G dG, then 

limg'^(r'^) = limqe('re) = q* ■ 

e -)-0 e ->-0 


(4.5) 


Proof. First, note that qe{t) is the first component of the 2(i-dimensional Markov process 
Ze{t) = {qe{t), pe{t)). Because of the structure of the p-component of the drift of this process 
and our assumptions on the vector field b, starting from {q,p) G the trajectory of Ze{t) 
spends most of the time in a small neighborhood of the point q = O and p = 0, with 
probability close tol, as0<e<<l. From time to time, the process Ze{t) deviates from 


this point and, as proven in Theorem 2.1, the deviations of qe{t) are governed by the large 


deviation principle with action functional /, defined in (2.4). This allows to prove the 


validity of (4.3), (4.4) and (4.5) in the same way as Theorems 4.41, 4.42 and 4.2.1 from m 
are proven. We omit the details. 

□ 


As an immediate consequence of (4.2) and m Theorem 4.3.1], we have the following 
result. 


Theorem 4.2. Assume a{q) := a{q)cr*{q) = I and a{q)b{q) = —VU{q) + l{q), for any 
q G for some smooth function U : —>• M having a unique critical point (a minimum) 
at O € and such that 

{VU{q),l{q))=Q, g G M'". 

Then 

V{q) = 2U{q), q G M-". 

From Theorems |4.1| and |4.2[ it is possible to get a number of results concerning the 
asymptotic behavior, as e | 0, of the solutions of the degenerate parabolic and the elliptic 
problems associated with the differential operator defined by 

1 ^ d‘^u ( 1 \ 

C^u{q,p) = 2^1 - -o:iq)pj ■S/pu{q,p) +p-S/gu{q,p). 


Assume now that the dynamical system (4.1) has several asymptotically stable attrac¬ 
tors. Assume, for the sake of brevity, that all attractors are just stable equilibriums Oi, 
O 2 ,... ,Oi. Denote by £ the set of separatrices separating the basins of these attractors, and 
assume the set £ to have dimension strictly less than d. Moreover, let each trajectory q{t), 
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starting at qo ^ \ S, be attracted to one of the stable equilibriums Oi, i = 1,... ,1, as 

t —>■ oo. Finally, assume that the projection of b{q) on the radius connecting the origin in 
and the point g G is directed to the origin and its length is bounded from below by 
some uniform constant 0 > 0 (this condition provides the positive recurrence of the process 
Zeit) = {qe{t),Peit)), t > 0). 

In what follows, we shall denote 

^(91,92) 


= - inf 


Oi{f{s))cr ^(/(s)) (j{s) - ds, /(O) = qi, f{T) = q2, T > O} 


and 


Vij = V{Oi,Oj), i,jG 

In a generic case, the behavior of the process (q^(t),p‘^(t)), on time intervals of order 
exp(Ae“^), A > 0 and 0 < e << 1, can be described by a hierarchy of cycles as in [TT] 
and [6]. The cycles are defined by the numbers Vij. For (almost) each initial point q and a 
time scale A, these numbers define also the metastable state Oi*, i* = i*{q, A), where q'^{t) 
spends most of the time during the time interval [0, exp(Ae“^)]. Slow changes of the field 
b{q) and/or of the damping coefficient a{q) can lead to stochastic resonance (compare with 

m)- 


Consider next the reaction diffusion equation in 

du, 


dt 


{t, q) = Cu{t, q) + c{q, u{t, q))u{t, q), 


(4.6) 


, u{0,q) = g{q), g G f > 0. 

Here £ is a linear second order uniformly elliptic operator, with regular enough coefficients. 
Let q{t) be the diffusion process in associated with the operator £. The Feynman-Kac 
formula says that u can be seen as the solution of the problem 

u{t,q) =Eqg{q{t))ex.-p f c{q{s),u{t - s,q{s)) ds. (4.7) 

40 

Reaction-diffusion equations describe the interaction between particle transport dehned 
by q{t) and reaction which consists of multiplication (if c{q, u) > 0) and annihilation (if 
c{q, u) < 0) of particles. In classical reaction-diffusion equations, the Langevin dynamics 
which describes a diffusion with inertia is replaced by its vanishing mass approximation. If 
the transport is described by the Langevin dynamics itself, equation ( |4.6| ) should be replaced 
by an equation in Assuming that the drift is equal to zero {b{q) = 0), and the damping 
is of order e~^, as e 0, this equation has the form 

du'' 1 1 

-^{t,q,p) = 9 X/ “*4(9)777:177:(9,P) - -a{q)p-Vpu''{q,p) + p ■Vqu‘'{q,p) 


dpidpj 


*,i=i 

+c{q,u^{t,q,p))u''{t,q,p), 
u^{0,q,p) = g{q) > 0, {q,p) G 


i > 0, (g,p) G 


h2d 


(4.8) 
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Now, we define 


where 


= sup c{f{s),0)ds - Itif) : /(O) = g, /(t) G Go 


= a^(/(s))a ^{f{s))f{s)- f{s)ds, 


and Go = supp{g(g), q G 

Definition 4.3. 1. We say that Condition (N) is satisfied if R{t,x) can be characterized, 

for any t > 0 and x £ T,t = {q £ R{t, q) = 0}, as 

sup c(/(s), 0) ds - Itif), /(O) = q, f{t) £ Go, R{t - s, f{s)) <0, 0 < s < 11 . 


2. We say that the non-linear term f{q,u) = c{q,u)u in equation (4.8) is of KPP 
(Kolmogorov-Petrovskii-Piskunov) type if c{q,u) is Lipschitz-continuous, c(g, 0) > 
c(g, u) > 0, for any 0 < u < 1, c{q, 1) = 0 and c{q, u) < 0, for any u > 1. 


Theorem 4.4. Let the non-linear term in (4.8) be of KPP type. Assume that Condition 
(N) is satisfied and assume that the closure of Gq = supp{( 7 (g), q £ M'^} coincides with the 
closure of the interior o/Go. Then, 


and 


limK{t/e,q,p) = 0, ifR{t,q)<0, 

e^O 


limK{t/e,q,p) = I, if R{t,q) > 0, 

e^O 

p2d 


(4.9) 


(4.10) 


so that equation R{t,q) = 0 in defines the interface separating the area where K, the 
solution of (4.8), is close to 1 and to 0, as e f 0. 


Proof. If we define ufit,q,p) = K{t/e,q,p), the analog of (4.7) yields 


Ue{t,q,p) = E(q_p)fi((g,(f))exp 


c{qfis),ufit - s,qfis),pfis)) ds 


(4.11) 


where zfit) = (qfis),pfis)) is the solution to equation (2.6). By taking into account our 
assumptions on c{q,u), we derive from (4.11) 


Ue{t,q,p) <E(^q^p)g{qfit))exp[- / c{qfis),0)ds 


Theorem |2 .1 1 and the Laplace formula imply that the right hand side of the above inequality 
is logarithmically equivalent , as e 0, to exp (^i?(t,g)) and this implies (4.9). 

In order to prove (4.10), first of all one should check that if R{t,q) = 0, then for each 
5 > 0 

Ue{t,q,p) > exp (--<5) , (4.12) 
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when e > 0 is small enough. This follows from (4.11) and Condition (N), if one takes into 
account the continuity of c{q,u). The strong Markov property of the process {q^{t), 
and bound (4.12) imply (4.10) (compare with [5]). 

□ 


Consider, as an example, the case c(g, 0) = c = const. Then 


R{t, q) =ct- inf {/*(/), /(O) = q, f{t) G Gq} . 


The infimum in the equality above coincides with 

^p^{q,Go), 


(see, for instance, [5] for a proof), where p{qi,q 2 ), qi,q 2 £ 1^^, is the distance in the 
Riemaniann metric 

d 


ds 



aij{q)dqidqj. 

ij=l 


This implies that the interface moves according to the Huygens principle with the constant 
speed if calculated in the Riemannian metric ds. 

If a(q) = 0 in a domain Gi C M'^, the points of Gi should be identified. The Riemaniann 
metric in induces now, in a natural way, a new metric p in this space with identified 
points. The motion of the interface, in this case, can be described by the Huygens principle 
with constant velocity in the metric p. 

If c{q,0) is not constant, the motion of the interface, in general, cannot be described 
by a Huygens principle. Actually, the motion can have jumps and other specific features 
(compare with [5]). 

Finally, if the Condition (N) is not satisfied, the function R{t, q) should be replaced by 
another one. Define 


R{t, q) = sup 


mm 

0<a<t 


c(/(s), 0) ds - laif) ■■ /(O) = q, f{t) £ Go 


The function R{t, q) is Lipschitz continuous and non-positive and if Condition (N) is satis¬ 
fied, then 

R{t, q) = min{i2(t, O'), 0} . 

By proceeding as in [6], it is possible to prove that 


lim M'^(t/e, g,p) = 0, ifi?(t,g)<0, 


and 


\miu\t/e,q,p) = 1, 


if (t, q) is in the interior of the set {{t,q) : t > 0, q £ R{t, q) = 0}. 


Finally, we would like to mention a few generalizations. 
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1. The arguments that we we have used in the proof of Theorem 2.1, can be used to 
prove the same result for the equation 


r (t) = ww) - + 4 

g"(0) = g e g"(0) = p € R‘^, 


for any f3 < 1/2. As a matter of fact, with the very same method we can show that 
also in this case the family {qe}e>o satisfies a large deviation principle in C'([0, T]; M'^) 
with action functional I and with normalizing factor 


2. The damping can be assumed to be anisotropic. This means that the coefficient a{q) 
can be replaced by a matrix a{q), with all eigenvalues having negative real part. 

3. Systems with strong non-linear damping can be considered. Namely, let {q^^p'") be the 
time-inhomogeneous Markov process corresponding to the following initial-boundary 
value problem for a degenerate quasi-linear equation on a bounded regular domain 
G C 


du^{t,p,q) 

dt 


I 

*, 1=1 
a{q,u'^{t,q,p)) 


q,p) 


dpidpj 


+ b{q) ■ Vpu\t,q,p) 


p-Vpu''{t,q,p) +p-Vqu''{t,q,p). 


. u"{0,q,p) = g{q), u^{t:Q,P)\^^aG 


Existence and uniqueness of such degenerate problem, under some mild conditions, 
follows from [U Chapter 5]. The non-linearity of the damping leads to some pecu- 
larities in the exit problem and in metastability. In particular, in the generic case, 
metastable distributions can be distributions among several asymptotic attractors and 
the limiting exit distributions may have a density (see HO]). 
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